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HELIUM IN THE TOPSIDE VENUS IONOSPHERE 


ABSTRACT 


The quantity of helium in the Venus atmosphere can be estimated 
from an examination of the measured ionization profiles. The amount 
of helium necessary to give agreement between the computed and exper- 
imental results is approximately 4 x 10 8 cm -3 at 140 km. When an 
eddy diffusion coefficient of 10 6 cm 2 sec -1 is used, the" atmospheric 
helium mixing ratio is found to be 6 x 10" 5 ; a value nearly ten times 
that for Earth. 


ir 


HELIUM IN THE TOPSIDE VENUS IONOSPHERE 

The best direct evidence for the presence of helium on Venus is ob- 
tained from the recent detection of its radioactive progenitors on the 
surface of Venus by Venera 8 (announced in Tass, 1972). Since the 
quantities of radioactive materials found were similar to measurements 
of Earth samples, the ultimate detection of helium on Venus should be a 
near certainty. Three years prior to the landing of Venera 8, Knudsen 
and Anderson (1969) assumed that the Venus helium production rate was 
similar to Earth's, but that the loss rate was due only to thermal es- 
cape (negligible relative to production), and calculated the helium atmos- 
pheric mixing ratio as 2 x 10 ^ Since this value is 40 times that for 
Earth, an attempt should be made to check for consistency with other data. 

An additional measure of the helium content is afforded by the upper 
atmosphere and ionosphere measurements obtained from Mariner V. The 
Lyman a determination of the hydrogen density in the exosphere, Barth (1968), 
and the electron and total neutral density from the radio occultation ex- 
periment, Fjeldbo and Eshelman (1969), can be used in conjunction with 
atmospheric modeling to place limits on the helium concentration. The 
principal features of these measurements are: (1) An upper limit on the 

hydrogen density of 5 x 10 3 cm -3 at 500 km altitude; (2) the observed ab- 
rupt change in the electron density scale height in the vicinity of 200 km; 

(3) the value of the electron density above 250 km (about 10 4 cm -3 ), and 

(4) the inferred neutral gas temperature in the isothermal region of about 


650 °K. 
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A reference neutral atmosphere for this work can be constructed by 
selecting boundary values at 100 km altitude for each neutral species 
that are consistent with the known measurements. To produce the desired 
calculated atmosphere, the energy, momentum, and continuity equations are 


integrated for 

both the neutral and ionic species. 

The neutral atmos- 

phere boundary conditions are adjusted to match the lower atmosphere 

measurements of 

total density and possible composition; approximately 

97% CO with up 

to 3% of minor constituents A r, N 

, 0 , H , He, CO, and 

2 

2 

2 2 

others, Marov (1972). Of these, the upper limits 

of N and 0 are about 

2 2 

2% and 0.1%, respectively. The remaining species 

are estimated at refer 

ence altitudes 

as shown in Table 1, and the complete atmospheric model 

in Figure 1. 




TABLE 1 




Reference 

SPECIES 

DENSITY (cm -3 ) 

Altitude (km) 

n 2 

5xlO n - 4xl0 12 

100 

C0 2 

2xl0 14 - 5x1 0 14 

100 

0 2 

lxlO 10 - 2xlO n 

100 

He 

0 - 2xl0 10 

100 

H 

<5x1 0 3 

500 

C 

0 - 2x1 0 7 

150 

0 

5x1 0 5 - 5x1 0 8 

150 

CO 

lxlO 9 - lxlO 11 

100 

N 

5x1 0 5 - 5x1 0 8 

150 

NO 

Unknown - <5x1 0 6 

100 



-3- 


The value of the total density at 100 km altitude was chosen by 
extrapolating the lower atmospheric data taken by the Mariner and Venera 
spacecraft. Ainsworth and Herman (1972), and by requiring that the 
electron density reaches a maximum near an altitude of 145 km for Mariner 
conditions. The base temperature at 100 km was varied between 160 and 
200°K, with 180°K adopted as a reference value. This temperature varia- 
tion is able to produce a shift in the peak altitude of the electron 
density of + 5 km. Its effect on the He and He + distribution is small, 
and therefore a base temperature of 180°K is used throughout. 

Since the purpose of this work is to discuss the effect of the varia- 
tion of neutral helium on the resulting electron density profile above 
200 km, the appropriate equations and chemical reactions are summarized 
below. A one dimensional diffusion model is adopted using the continuity 
and momentum equations 


an 

_L + L (n v ) = P - L n (1) 

3z i i i i i 


av ap 

p i_ = - i_ - p g + en E - E K (v -v ) - z K v (2) 

i at az i i j ij i j j ij i 

IONS NEUTRALS 

av 

where the higher order terms p v i_ , thermal diffusion, and the motion 

i i az 

of the neutral atmosphere are neglected. The possibly unfamiliar symbols 

are K = n y v , where y is the reduced mass and v is the collision 
ij i ij ij ij ij 

frequency, Banks (1966), Dalgarno (1969), P is the total volume production 

i 
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rate from both solar and chemical sources for the i th ionic species, 

and L n is the corresponding chemical loss rate. In the case of helium, 
i i 

the total production of He is from the solar EUV ionization of He com- 
puted in the usual manner, Herman et al . (1971), based on the solar flux 
given by Hinteregger (1970), and that arising from energetic secondary 
photoelectrons. The chemical loss channels are described in Table 2. 

The loss channels via atomic species all appear to have very small loss 
rates, Fehsenfeld et al. (1966). It is that fact and the low measured 
value of [H] that is mainly responsible for He + being the dominant top- 
side ion. Much of the H + produced is lost by charge exchange with 0 at 

3. 0 3 ! 

a rate of 3.75 x 10 cm" sec” as well as by other rapid losses to 
CQ 2 and 0 2 » Banks (1972). Another possibil tty. for a topside Venus ion 
is C + since it is produced rapidly both by photoionization and charge 
transfer processes, and is not lost at a large rate to other atomic 
species. However, unless there are very large amounts of neutral carbon 
present, its relatively heavy mass compared to He precludes it from a 
significant contribution to the topside ions. Most of the C + that is pro- 
duced is lost to CO 2 and O 2 with rates of 3.7 x 10" 9 and 2.2 x 10" 9 cm" 3 sec -1 
respectively, McDaniel et al . (1970). 

The result of using the above reaction rates, the appropriate photo- 
ionization cross sections, Hinteregger (1965), McDaniel (1964), and 
Sullivan and Holland (1964), and. equations (1) and (2), yield the results 
shown in Figure 2. The major ions adding up to the electron density 

n are indicated in four altitude regimes. The bottomside is composed 
e 

mainly of C02 + , 0 2 + , and N0 + . The relative amounts are governed 
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TABLE 2 


Reaction 

Rate 

Reference 

H + +0 o -*H +0+0 + 
e L e 

1.5(-9> 

Fehsenfeld (1966) 

H e ++ V°2 ++H e 

3 . 0( -10) 

Warneck (1966) 

H e + +C0 2 +C0 + +0+H e 

9.6(-10) 

Fehsenfeld (1966) 

-K) + +C0+H e 

1.8(-10) 

Fehsenfeld (1966) 

->C0 2 + +H e 

6.0(-l 1 ) 

Fehsenfeld (1966) 

H + +CCH-H +0+C + 
e e 

1.7(-9) 

Fehsenfeld (1966) 

H + +N?-*-H +N+N + 
e ^ e 

7.1 (-10) 

Warneck (1966) 

* H e +N 2 + 

7 . 6 (-1 0) 

Warneck (1966) 

H + +e->-He+hv 
e 

5 . 0 ( -1 2 ) 

Biondi (1964) 


by the amount of N and 0 present in the atmosphere. N0 + is produced 

copiously, even in the absence of NO, and can amount to half the total 

ion density at the electron density maximum based on the limits of the 

neutral atmosphere given in Table 1. As both N and 0 decrease with 

altitude the amount of N0 + also decreases rapidly giving way to 0 , CO^, 

and 0 + . With further increases in altitude He + dominates all other 
2 

ions present. 

The observations of n from Mariner V show a topside ion density 

e 

of about 10 + cm" 3 and a very abrupt slope change between the molecular 
ion to the atomic ion region. Comparing the topside data to the results 
in Figure 2 show that the measured profile falls between cases (b) and (c). 
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Accordingly, limits can be placed on the helium density above the tur- 

8 8 

bopause of from 3.5 x 10 to 5 x 10 at 140 km altitude. When an 

6 2 f 

eddy diffusion coefficient of K = 10 cm sec" is used, the helium 

10 10 _3 

density well below the turbopause if from 1x10 to 1.5 x 10 cm at 
100 km and leads to cases (b) and (c). This corresponds to a mixing 

-5 

ratio of approximately 6 x 10 . The results obtained here, which in- 

clude the effects of ion diffusion and self consistency with the tem- 
perature calculation, require significantly more neutral helium than 
previous estimates. Whitten (1970) obtained an estimate of about 0.1 
of the present value from a dayside ionospheric calculation. McElroy 
and Strobel (1969) estimated about 0.1 from a nightside calculation and 
Knudsen and Anderson (1969) about 1/4. These factors are significant 
because a reduction by a factor of 1/2 produces the results shown in 
case (a) with a drastically decreased scale height and too small a value 
of electron density. In the case of Whitten, enough approximations 
were made in evaluating the terms in the momentum and continuity equations 
that the accuracy of his results are difficult to evaluate. In particu- 
lar, one of his key results is that the scale heights of C0 + and 0 + are 
the same as He + . The results obtained from solving for each species shows 

+ 3 _3 

that CO reaches a maximum of about 10 cm at 250 km altitude where 
the dominant ions are 0 2 + , C0 2 + , 0 + , and He + and decreases rapidly above 

250 km, (C0 + = 40 at 350 km). Similarly, 0 + reaches a maximum of about 

3 3 

8 x 10 at 260 km and decreases to about 1 xlO at 360 km. The very 

different He + scale height can be seen in Figure 2. 
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-4 

McElroy and Strobe! (1969) estimated a mixing ratio of 10 by 
making an Earth analogy for the turbopause height. Since a lowering of 
the turbopause altitude (a strong function of total density) results in 
a lower requirement for the mixing ratio, the value estimates here, 

_ 5 

6 x 10 , does not contradict their method of estimation. 

The method employed here for determining where the cessation of 
turbulent mixing and the onset of diffusive separation occurs is based 
on the relative values of K and D, the eddy and binary diffusion coeffic- 

o 

ients, respectively. For altitudes where K >D, K = K , and in the altitude 

o o 

range where K ^ D or K <D,K decreases rapidly so that the onset of 

o o 

diffusive separation takes place within about a 5 to 10 km range. This 
is verified experimentally for Earth, Keneshea and Zimmerman (1970) 
with the results approximated quite well by the form K * K exp [-D/K ]. 

o o 

Expressing K in this manner permits their result to be scaled to the 
Venus atmosphere. 

If, as McElroy and Strobe! suggest, the nightside ionization above 
200 km is to be maintained by lateral flow from the dayside, the amount 
of helium in the present model is sufficient. It should be noted that 

4 3 

the nightside ionization ledge of 10 cm below 200 km cannot be main- 
tained by flow from the dayside because of the rapid He + + C0 2 charge 
exchange rate. In addition, if such a process were considered, it would 
require downward diffusion from higher altitudes. This He + flow in the 
presence of C0 2 + ions would not be possible due to the presence of an 
upward electric field resulting from a C0 2 + - He + mixture. Perhaps a 
combination of their suggestions, with the ionization ledge below 
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200 km produced by an influx of energetic solar wind particles, and 
with the topside produced by He + flow from the dayside would be appropri 

ate. ' ' • v 

In summary, the observed dayside ionization above 200 km csn.be 

satisfactorily accounted for with a helium density of about 4 x 10 8 cm -3 
at 140 km altitude. This corresponds to a mixing ratio of 5 x 10 if 
the eddy diffusion coefficient is taken to be 10 cm sec * Such a 
model is consistent with estimates of helium generated from radioactive 
sources over geological time, and with the requirements for maintenance 
of the nightside ionosphere above 200 km by lateral flow of He from 
the dayside of Venus. 
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CAPTIONS 


1: A model of the Venus neutral atmosphere above 100 km altitude. 

Dashed lines refer to the upper scale and solid lines to 
the lower scale. Thick lines are molecular species and thin 
lines are atomic species. The boundary values are based on 
Table 1, and for all but CO 2 and H are upper limits esti- 
mated from lower atmospheric measurements or from chemical 
equilibrium. 


2: The calculated Venus ionosphere showing the major ions in 

their respective regions of dominance. The topside iono- 
sphere is shown for four cases of helium density based on 
boundary values at 100 km. Solid lines are the electron 
density and the dashed lines are the H e density. Two 
cases are shown for the bottomside electron density for two 
different neutral gas temperatures T = 160°K and T = 180°K 
at 100 km. 
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